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We examined morphological differentiation and allozyme diversity in nine populations of Vincetoxicum 
katoi and its close relative V yamanakae. Although, V yamanakae has been distinguished from V. katoi 
in having a larger corolla and gynostegium, the principal component analysis based on 12 floral and three 
leaf characters did not support the distinction of the two taxa. The six allozyme loci examined revealed 
that population of V katoi and V yamanakae maintained high genetic diversity (P = 72.1, A = 2.40, h = 
0.313). The nine populations were clustered into three regions (Kanto-Tokai; middle Kiniki; and Kii- 
Shikoku) based on genetic differences, contradicting the current taxonomic treatment. Considering the 
evidence, it is appropriate to treat V. vamanakae as a synonym of V. katoi. The high genetic differentia¬ 
tion among regions suggested that disjunct distribution of the V. katoi-V. yamanakae complex might re¬ 
flect the persistence of refugia since the last glacial period. In particular, the middle Kinki (Hyogo) pop¬ 
ulation is located outside of known evergreen forest refugia, suggesting that it might have survived dur¬ 
ing the glacial period within cooler vegetative zones, such as in temperate forests, whereas the popula¬ 
tions on the Pacific-side retreated to warm-temperature coastal refugia. 

Key words: allozyme, disjunct distribution, floral morophology, Vincetoxicum katoi , 

Vincetoxicum yamanakae 


The Japanese flora is affected not only by the 
current climate but also by past geographic events 
(Aoki et al. 2006). Recent phylogeographic stud¬ 
ies of plants in alpine-sub alpine regions (e.g., Fu- 
jii & Senni 2006, Ikeda et al. 2008), temperate 
forests (e.g., Hiraoka & Tomaru 2009, Kikuchi et 
al. 2010), and evergreen forests (e.g., Aoki et al. 
2004, 2006) have provided new insights into the 
location of glacial refugia and the routes of post¬ 
glacial range expansion of plants within the Japa¬ 
nese Archipelago. Palynological records for low¬ 


land evergreen forests in Japan have suggested 
the existence of multiple refugia along the Pacific 
coast during the last glacial period (Kamei and 
Research Group for Biogeography from Wiirrn 
Glacial 1981, Hattori 2002). Using chloroplast 
DNA sequence variations, Aoki et al. (2004, 
2006) confirmed the existence of refugia for sev¬ 
eral evergreen tree species outside of Kyusyu. 
The locations of some refugia, however, were 
only weakly confirmed by molecular-based anal¬ 
yses. 
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Disjunct distributions, the presence of plants 
or animals occurring in two or more remote geo¬ 
graphic areas, are thought to be caused by three 
major processes (Raven 1972, Cox & Moore 
2005, Lomolino et al. 2005): 1) ancient tectonic 
events such as continental drift that resulted in 
intercontinental splits in species distributions, 2) 
more recent climate changes (the glacial-intergla¬ 
cial cycle during the Quaternary), and 3) long¬ 
distance dispersal. Alpine plants, which exhibit 
island-like disjunct distributions, have been the 
subject of intensive research on glacial refugia 
using molecular phylogenetic approaches (e.g., 
Fujii & Senni 2006). As in alpine regions, low¬ 
land plants with disjunct distributions are expect¬ 
ed to provide information regarding drastic dis¬ 
tribution changes during glacial periods. This is 


because even though lowlands have been affected 
by long-term deforestation, their current disjunct 
distributions may represent footprints of glacial 
refugia. 

Vincetoxicum katoi (Ohwi) Kitagawa is an en¬ 
demic perennial understory herb that occurs 
chiefly in Cryptomeria japonica (L. fil.) D. Don 
plantations, and is disjunctly distributed from 
lowlands to hills on the Pacific side of the Kanto 
district and in the southern and middle parts of 
the Kinki (Kii peninsula and Hyogo) district on 
Honshu, Japan (Fig. 1) (Yamazaki 1993). Vince¬ 
toxicum yamanakae (Ohashi & Ohwi) Ohashi is 
closely related to V katoi, and also occurs chiefly 
in Cryptomeria japonica plantations in the low¬ 
lands of Shikoku (Fig. 1) (Yamazaki 1993). V ya¬ 
manakae is distinguished from V. katoi by having 


K3 K1 



Fig. 1. Geographic distribution of the V. katoi (open circle)-K yamanakae (closed circle) complex based on herbarium speci¬ 
men information and location of sampling sites. For population codes, see Table 1. 
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a larger corolla and gynostegium (Yamazaki 
1993). Currently, V. katoi is listed as ‘vulnerable’, 
whereas V vamanakae is listed as ‘critically en¬ 
dangered’ because of its more restricted distribu¬ 
tion (Environment Agency of Japan 2007). As 
suggested by Yamazaki (1993), however, the co¬ 
rolla size in species of asclepiads is variable. Am¬ 
biguity in taxonomic status may lead to errone¬ 
ous conservation decisions and designations 
(Frankham et al. 2004). To assign the proper con¬ 
servation designation to species such as V. katoi 
and V vamanakae , their morphological and ge¬ 
netic distinction should be determined. 

The first objective of this study was to ascer¬ 
tain the distinction between Vincetoxicum katoi 
and V yamanakae by quantitatively evaluating 
their morphological traits. The second objective 
was to assess the genetic diversity and differen¬ 
tiation of populations of V. katoi and V. ya¬ 
manakae using allozymes. We also discuss the 
factors that influence the current disjunct distri¬ 
bution of the populations of V katoi and V ya¬ 
manakae. 

Materials and Methods 

Plant materials 

Nine populations of Vincetoxicum katoi and 
V. yamanakae (Fig. 1, Table 1) were sampled and 
tentatively classified as either V. katoi (designated 
K) or V yamanakae (designated Y) following 
Yamazaki’s (1993) report of their geographical 
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distributions. Five populations from Kanto, To- 
kai, Kii, and middle Kinki were classified as V 
katoi, and four populations from Shikoku were 
classified as V. yamanakae. A combination of in¬ 
florescences and leaves was collected from flow¬ 
ering individuals. Inflorescences were fixed in 
FAA (a mixture of formalin, glacial acetic acid, 
and 70% ethanol at a ratio of 5:5:90 by volume) 
for measurement of floral traits. Leaves were 
transported on ice to the laboratory and refriger¬ 
ated (4°C) for several days until being extracted 
for electrophoresis analysis. Voucher specimens 
representing each population were deposited in 
the Flerbarium of Tohoku University (TUS), 
Shoei Junior College (SHO) and Tokushima Pre- 
fectural Museum (TKPM). 

Morphology 

Morphological measurements were recorded 
for 72 and 90 individuals from four populations 
each of Vincetoxicum katoi and V. yamanakae, 
respectively (Table 1). A digital caliper and a mi¬ 
crometer under a dissecting microscope were 
used to measure 12 floral features: pedicel length 
(PL), corolla diameter (CLD), corolla lobe length 
and width (CLL and CLW, respectively), corona 
diameter (CND), length of corona lobe (CNL), di¬ 
ameter and length of the gynostegium (GYD and 
GYL, respectively), length and width of the cor- 
pusculum (COL and COW, respectively), and 
length and width of the pollinia (POL and POW, 
respectively) (Fig. 2). For leaf samples, three 


TABLE 1. Population code, localities, altitude and sample size for morphological (N M ) and allozyme (VO analysis for the popu¬ 
lations examined in the V. katoi-V. yamanakae complex. 


Population 

code 

Locality 

Altitude 

n m 

/v A 

Voucher* 

K1 

Chiba Pref., Tomisato 

40 m 

25 

24 

T.Y.&A.Y. 7416 { TUS) 

K2 

Shizuoka Pref.. Fukuroi 

40 m 

- 

13 

T.Y.&A.Y. 7181 (TUS) 

K3 

Aichi Pref., Shinshiro 

200 m 

19 

31 

T.Y.&A.Y. 7421 (TUS) 

K4 

Wakayama Pref., Shingu 

60 m 

8 

6 

T.Y.&A.Y. 11169 (TKPM) 

K5 

Hyogo Pref., Kobe 

600 m 

20 

25 

T. K. 28983 (SHO) 

Y1 

Tokushima Pref., Kaiyo 

60 m 

24 

24 

T. Y. 11175 ( TKPM) 

Y2 

Tokushima Pref., Kaiyo 

70 m 

23 

27 

T. Y. 11183 (TKPM) 

Y3 

Kochi Pref., Konan 

40 m 

25 

51 

T.Y.&A.Y. 7438 (TUS) 

Y4 

Kochi Pref., Tosashimizu 

140 m 

18 

25 

T.Y.&A.Y. 7198 (TUS) 


* T. Y. ; Yamashiro, T.; A. Y. ; Yamashiro, A.; T. K.; Kobayashi, T. 
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CLD 



Fig. 2. Floral characters measured in this study. A: whole flower, B: corona and gynostegium, C: pollinaria, PL: pedicel length, 
CLD: corolla diameter, CLL: corolla lobe length, CLW: corolla lobe width, CND: corona diameter, CNL: corona lobe 
length, GYL: gynostegium length, GYD: gynostegium diameter, COL: corpusculum length, COW: corpusculum width, 
POL: pollinia length, POW: pollinia width. 


characteristics, leaf blade length (LBL) and width 
(LBW), and petiole length (PET), were measured. 

Principal component analysis (PCA) was per¬ 
formed on 12 floral and three leaf characteristics 
to summarize the overall tendencies of the mor¬ 
phological variation using Excel-Tokei 2000 (So¬ 
cial Survey Research Information Co., Ltd., To¬ 


kyo, Japan). 

Allozyme electrophoresis 

For each plant, 100 mg of fresh leaf tissue was 
homogenized in lmL of an extraction buffer 
comprising 93 mM Tris-HCl (pH 7.5), 23.4% 
glycerol, 0.5% (v/v) Tween 80, 2.8 mM dithioth- 
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reitol, and 0.5% 2-mercaptoethanol (Uchida et al. 
1991). The homogenates were centrifuged at 
20,000 x g at 4°C for 15 min. The resulting super¬ 
natants were used as crude enzyme extracts. 
Polyacrylamide vertical slab gel electrophoresis 
was carried out according to procedures de¬ 
scribed by Davis (1964) and Ornstein (1964). 10 
juL crude enzyme extracts were analyzed by 
polyacrylamide gel electrophoresis. Six enzyme 
systems were examined: aspartate amino-trans¬ 
ferase (AAT, EC 2.6.1.1), colorimetric esterase 
(EST, EC 3.1.1), leucine aminopeptidase (LAP, 
EC 3.4.11.1), malate dehydrogenase (MDH, EC 
1.1.1.40), phosphoglucoisomerase (PGI, EC 
5.3.1.9), and shikimate dehydrogenase (SKDH, 
EC 1.1.1.25). The staining protocols used were 
those of Tsumura et al. (1990). 

The following indices were used to quantify 
genetic diversity within each population: the pro¬ 
portion of polymorphic loci ( P) at the 95% crite¬ 
rion, the number of alleles per locus (A), the ex¬ 
pected heterozygosity (/?), and the inbreeding co¬ 
efficient (F| S ). As in Hamrick and Godt (1989), 
the loci that were polymorphic in at least one 
population were treated as polymorphic at the 
species level. The significance of deviation of F IS 
from zero was evaluated using a permutation test. 
To calculate A, h , and F 1S , we used FSTAT 2.9.3 
(Goudet 2001). 

The genetic differentiation among the popula¬ 
tions of Vincetoxicum katoi and V yamanakae 
was determined using F ST of Weir & Cockerham 
(1984). The significance of pair-wise F ST values 
was tested by applying the adjusted sequential 
Bonferroni correction based on 1,000 permuta¬ 
tions. All estimates of F SJ and their tests of sig¬ 
nificance were performed using FSTAT 2.9.3 
(Goudet 2001). To access the overall population 
genetic differentiation tendency, the gene diver¬ 
sity statistics, corrected for small sample size er¬ 
rors, were also estimated (Nei & Chesser 1983). 
The genetic relationships among the V. katoi and 
V. yamanakae populations were evaluated by 
generating a neighbor-joining tree based on the 
unbiased genetic distance (Nei 1978) using POP- 
TREE2 (Takezaki et al. 2010). The significance 
of the topology was estimated with 1,000 boot¬ 


strap replicates. To assess components of genetic 
variation within and between geographical re¬ 
gions, analyses of molecular variance (AMOVAs; 
Excoffier et al. 1992) with a significances test of 
999 random permutations were performed using 
GenAlEx 6 (Peakall & Smouse 2006). 

Results 

Morphology 

The morphological characteristics of the eight 
populations analyzed are summarized in Appen¬ 
dix 1. The variances accounted for by PCI and 
PC2 were 46.2% and 11.0% of the morphometric 
data, respectively. The principal component load¬ 
ings on the first and second components are listed 
in Table 2. The plots of samples on PCI and PC2 
are shown in Fig. 3. The PCI represents well the 
total size of flowers and leaves. For the PC2, 
whereas the COL, COW, PL and POW were posi¬ 
tively correlated, the CNL, LBL, LBW and PET 
were negatively correlated. Therefore, samples 
having large pollinia, longer pedicels, smaller 
leaves and smaller corona tended to appear in the 
positive field for PC2. The K3 population clus¬ 
tered separately, although other populations were 
not clearly separated. 


Table 2. Principal component loading on the first and second 


components. 


Character 

Principal component 

1 

2 

PL 

0.400 

0.623 

CLD 

0.866 

0.037 

CLL 

0.802 

-0.098 

CLW 

0.793 

-0.068 

CND 

0.887 

-0.189 

CNL 

0.782 

-0.248 

GYD 

0.877 

0.151 

GYL 

0.889 

-0.090 

COL 

0.560 

0.368 

COW 

0.406 

0.627 

POL 

0.726 

0.052 

POW 

0.512 

0.241 

LBL 

0.536 

-0.352 

LBW 

0.324 

-0.225 

PET 

0.349 

-0.580 
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Fig. 3. Plot of principal component analysis for nine populations in the V. kaioi V. yamanakae complex. 


Genetic diversity 

A total of 10 putative loci were detected: Aat- 

1, Aat-2 , Est-1, Est-2, Lap, Mdh-1, Mdh-2, Pgi-1, 
Pgi-2, and Skdh. However, because the Aat-2 , Est- • 

2, Mdh-2, and Pgi-2 loci were not interpretable in 
all nine populations examined due to weak stain¬ 
ing, they were discarded. Among the six loci 
scored, Mdh-1 was monomorphic in all popula- 


Table 3. Proportion of polymorphic loci (/ > ), mean number 
of alleles (A), expected heterozygosity (/;), and inbreed¬ 
ing coefficient (F , s ) at six loci for the populations exam¬ 
ined. All F| S values did not significantly deviate from 
zero by permutation test. 


Population 

P 

A 

h 

As 

Kl 

83.3 

3.67 

0.488 

0.181 

K2 

83.3 

2.67 

0.443 

-0.005 

K3 

83.3 

3.00 

0.356 

0.022 

K4 

50.0 

1.67 

0.233 

-0.071 

K5 

66.7 

2.00 

0.169 

0.098 

Yl 

83.3 

2.17 

0.283 

0.021 

Y2 

50.0 

1.67 

0.174 

0.009 

Y3 

83.3 

2.83 

0.373 

0.118 

Y4 

66.7 

2.00 

0.300 

-0.034 

Mean 

72,1 

2.40 

0.313 

0.037 

Overall 

83.3 

4.17 

0.468 

0.119 


tions. The other five loci were polymorphic in at 
least one population. Allele frequencies of the six 
loci examined for each population are shown in 
Appendix 2. Table 3 shows the P, A, h and F IS val¬ 
ues for each population and for the overall popu¬ 
lation. The percentage of polymorphic loci (P) 
ranged from 50.0% (K4 and Y2) to 83.3% (Kl-3, 
Yl, and Y3). The mean number of alleles (A) 
ranged from 1.67 (K4 and Y2) to 3.67 (Kl). The 
expected heterozygosity (h) ranged from 0.169 
(K5) to 0.488 (Kl). The genetic diversity was 
highest in the Kl population and lowest in the Y2 
population. In none of the nine populations ex¬ 
amined did the inbreeding coefficient (F IS ) value 
deviate significantly from zero (Table 3). 

Genetic structure 

Pairwise comparisons of genetic differentia¬ 
tion (F ST ) across all six loci in the nine popula¬ 
tions examined are shown in Table 4. The F ST val¬ 
ues between the populations ranged from 0.042 to 
0.766, with an average of 0.337. The K5 popula¬ 
tion was highly genetically differentiated from 
the other populations (0.395-0.766). The genetic 
differentiation between the Kanto-Tokai and 


Nil-Electronic Library Service 




The Japanese Society for Plant Systematics 


October 2012 Yamashiro & al. —Population Genetic Diversity in Vincetoxicum 35 


Table 4. Pair-wise F ST among nine populations of V katoi-V. vamanakae complex. All values significantly departed from zero. 


Population 

K2 

K3 

K4 

K5 

Y1 

Y2 

Y3 

Y4 

Kl 

0.090 

0.152 

0.271 

0.395 

0.323 

0.468 

0.284 

0.327 

K2 


0.042 

0.271 

0.451 

0.292 

0.529 

0.291 

0.280 

K3 



0.322 

0.464 

0.343 

0.525 

0.308 

0.273 

K4 




0.640 

0.085" 

0.445 

0.082" 

0.175 

K5 





0.597 

0.766 

0.495 

0.596 

Y1 






0.336 

0.130 

0.215 

Y2 







0.243 

0.428 

Y3 








0.206 


#indicates non significant after sequential Bonferroni correction. 


Kii-Shikoku regions was also large (F ST = 0.271— 
0.529). The genetic differentiation among the 
Kii-Shikoku populations was higher than in the 
Kanto-Tokai populations, with the exception of 
K4-Y1, and K4-Y3 (not significant after Bonffer- 
oni correction). The total genetic diversity (// T ) 
and coefficient of gene differentiation (0 ST ) for 
all populations of Vincetoxicum katoi and V ya- 
manakae was 0.488 and 0.356, respectively. The 
latter value indicates that 64.4% of the total ge¬ 
netic variation was shared by all nine populations 
and the remaining 35.6% was partitioned within 


each population. 

A neighbor-joining tree based on Nei’s (1978) 
unbiased genetic distance for each population is 
shown in Fig. 4. The nine populations examined 
were largely divided into two clusters. One in¬ 
cluded the Kl, 2, 3, and 5 populations, and the 
other comprised the Yl-4 and K4 populations. 

Due to the NJ population clustering result, we 
performed an AMOVA analysis of the nine popu¬ 
lations after division into three regions (Kl-3, 
K5, Yl-3 with K4). The genetic divergence be¬ 
tween the three regions was highly significant 



Y2 


K5 


Fig. 4. Phenogram using neighbor-joining method based on unbiased genetic distance (Nei 1978) for nine populations in the V 
katoi-V. yamanakae complex. 
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Table 5. Results of analysis of molecular variance (AMOVA) for populations of the V katoi-V. yamanakae complex in three 
regions (Kanto-Tokai, middle Kinki, and Kii-Shikoku). 


Source of 

df 

Sum of 

variance 

% of 

Phi 

P 

variation 


squares 

components 

variation 



Among regions 

1 

284.432 

1.825 

36% 

PhiRT = 0.364 

<0.001 

Among populations 

7 

151.340 

0.969 

19% 

Phi PR = 0.304 

<0.001 

within regions 

Within populatios 

218 

483.778 

2.219 

44% 

PhiPT= 0.557 

<0.001 

Total 

226 

919.551 

5.013 

100% 




(PhiRT= 0.282, P < 0.001) (Table 5). 

Discussion 

Taxonomic status ofV. katoi and V. yamanakae 

Vincetoxicum katoi had been previously dis¬ 
tinguished from V yamanakae by its smaller co¬ 
rolla and corona (Yamazaki 1993). In contrast, 
our quantitative analysis of floral morphology did 
not support the discontinuity of these species. 
Only the K3 population was separated from the 
other populations. Certainly, some herbarium 
specimens collected in the Tokai region exhibit 
small flowers (corolla lobes 2.5 4 mm in length) 
and leaves in good agreement with the original 
description of V katoi (Ohwi 1943). Other popu¬ 
lations in the V. katoi distribution area, however, 
were not morphologically distinguishable from V 
yamanakae. Furthermore, because the Kii penin¬ 
sula population clustered with populations from 
Shikoku, the neighbor-joining tree based on six 
allozyme loci did not support the recognition of 
two species. Thus, it is appropriate to treat V. ya¬ 
manakae as a synonym of V. katoi. 

Genetic diversity 

Despite their highly restricted and disjunct 
distribution, our allozyme analysis revealed that 
members of the V. katoi-V. yamanakae complex 
maintain high genetic diversity at both the spe¬ 
cies and population levels. At the population lev¬ 
el, the three genetic diversity parameters of V. ka¬ 
toi-V. yamanakae complex (P = 72.1, A = 2.40, h 
= 0.313) were higher than those of endemic spe¬ 
cies {P = 41.3, A = 1.80, h = 0.069) and perennial 
herbs (P = 41.3, A = 1.79, h = 0.012) (Hamrick & 


Godt 1989). When compared with a highly out- 
crossing North American species of Asclepias (A. 
exaltata L.: A = 2.36, h = 0.182, Broyles & Wyatt 
(1993)), members of the V. katoi-V. yamanakae 
complex exhibit higher genetic diversity. In the 
nine populations examined, inbreeding coeffi¬ 
cient values did not significantly deviate from 
zero, indicating that outcrossing predominates 
(Table 1). Outcrossing species tend to retain 
greater genetic variation than selfing species 
(Loveless & Hamrick 1984). Therefore, outcross¬ 
ing may contribute to the maintenance of high ge¬ 
netic diversity in the V katoi-V. yamanakae com¬ 
plex. Another possible factor in the high genetic 
diversity is species of ancient origin and with a 
long evolutionary history (e.g., Kirengeshoma: 
Zhang et al. 2006, Chang et al. 2007). Diversifi¬ 
cation in the Japanese endemic species of Vince¬ 
toxicum , however, is thought to have occurred 
relatively recently due to the paucity of cpDNA 
and nrDNA variability (Yamashiro et al. 2004, 
2008). The high genetic diversity in V. katoi-V. 
yamanakae complex may be the result of ances¬ 
tral polymorphisms. Because the habitat in which 
they now occurred may have been altered to 
Cryptomeria japonica plantations more or less 50 
years ago, they might still maintain high genetic 
diversity. 

Genetic differentiation 

In general, outcrossing species are expected 
to exhibit limited genetic differentiation among 
populations due to increased gene flow (Loveless 
& Hamrick 1984). The coefficient of genetic dif¬ 
ferentiation (G st ) among the nine populations of 
the V. katoi—V. yamanakae complex was 0.356, 
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much higher than the mean for endemic plants 
(G st = 0.248), short-lived perennial herbs (G ST = 
0.233), and outcrossing animal-pollinated plants 
(<7 st = 0.197) (Hamrick & Godt 1989). The ge¬ 
netic differentiation among the populations of the 
V katoi-V. yamanakae complex greatly exceeded 
the genetic differentiation in the species of Ascle- 
pias in North America (Broyles & Wyatt 1993, 
Edwards & Wyatt 1994, Broyles 1998). Species of 
Asclepias are mainly pollinated by bumble-bees, 
leading to pollen-mediated long-distance gene 
flow, and therefore low population genetic differ¬ 
entiation (Broyles & Wyatt 1993). In contrast, 
pollination by fungus gnats was observed in the 
Kanto and Tokai populations of V katoi-V. ya¬ 
manakae complex (Yamashiro et al. 2008), sug¬ 
gesting that pollen-mediated gene flow is restrict¬ 
ed in the V katoi-V. yamanakae complex. Al¬ 
though these species have seeds with hairs (coma) 
that are likely wind-dispersed, they predominant¬ 
ly inhabit the closed understory of Cryptomeria 
japonica plantations, making their dispersal effi¬ 
ciency lower than in congeneric open-habitat spe¬ 
cies. In fact, the distribution of plants in the V 
katoi-V yamanakae complex is limited even on 
the same mountain; no other populations were 
found in neighboring valleys or in neighboring 
Cryptomeria japonica plantations. Colonization 
by seed dispersal may also be limited in these 
plants. 

The nine populations in the V katoi-V. ya¬ 
manakae complex were genetically divisible into 
three regions: Kanto-Tokai, middle Kinki (Hyo- 
go), and Kii-Shikoku. The genetic differentiation 
among these regions was relatively large ( PhiRT 
= 0.364 in AMOVA). 

The current regional genetic differentiation 
and distribution of plants in the V katoi-V. ya¬ 
manakae complex may be due to historical fac¬ 
tors. At the last glacial maximum (25,000-15,000 
years ago), Honshu and Shikoku were covered 
with boreal coniferous forests. Evergreen plants 
in Japan are presumed to have reestablished their 
postglacial range expansion from southern coast¬ 
al refugia during the last 12,000 years (Tsukada 
1984). The Pacific coast is thought to have been a 
refugium for evergreen forests during the most 


recent glacial maximum (Hattori 2002). In V ka¬ 
toi-V. yamanakae complex, the Kanto-Tokai and 
Kii-Shikoku, regions are located in close prox¬ 
imity and appear to have been glacial refugia, 
based on palynological analysis (e.g., Muroto, 
Ashizuri, Shionomisaki, Boso; Hattori 2002). 
The population differentiation within the Kanto- 
Tokai region (F ST = 0.042-0.152) is smaller than 
in the Kii-Shikoku region (0.130-0.428), sug¬ 
gesting that the Kanto Tokai populations may 
have expanded from a single refugium, whereas 
the Kii-Shikoku populations were derived from 
multiple refugia. 

In particular, the K5 population is located 
outside of known evergreen forest refugia. The 
post-glacial migration route on the Setouchi Sea 
side of the Kinki district was considered a route 
from Shikoku or the Kii peninsula (Hattori 2002, 
Aoki et al. 2004). Although the K5 population is 
genetically similar to populations in the Kanto- 
Tokai region than those in the Shikoku-Kii re¬ 
gions, long-distance dispersal of the Kanto-To¬ 
kai population to middle Kinki (Hyogo) is unlike¬ 
ly to have occurred, based on the large genetic 
differences between these regions (F ST = 0.395- 
0.464). Such a large genetic difference may indi¬ 
cate the existence of refugia for the V. katoi- V 
yamanakae complex around the K5 population, 
although evergreen forests are unlikely to have 
existed nearby. Our results suggest that not all 
species currently growing in the evergreen forest 
zone may have migrated in accompaniment with 
evergreen forest species. The K5 population may 
have survived during the glacial period within a 
cooler vegetative zone, such as in a temperate for¬ 
est, while the Pacific-side populations survived in 
coastal warm-temperature refugia. 

Taxonomy 

Vincetoxicum katoi (Ohwi) Kitag. J. Jap. Bot. 
34: 364 (1959). -Cynanchum katoi Ohwi, Acta 
Phytotax. Geobot. 12: 112 (1943). 

Cvnanchum yamanakae Ohwi & Ohashi, J. Jap. 
Bot. 48: 371 (1973), syn. nov. -Vincetoxicum ya¬ 
manakae (Ohwi & Ohashi) Ohashi, J. Jap. Bot. 
65:277 (1990). 
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Specimens examined'. JAPAN. Honshu. Chiba Pref.: Yot- 
sukaido, 19 May 1983, J. Haginiwa (OSA); Kamogawa, 
1933, T. Makino (MAK). —Shizuoka Pref.: Fujinomiya, 5 
Jun. 1940, J. Sngimoto (KYO); Kikugawa, 29 May 1932, 
G. Hashimoto (TI); Kanaya, 8 Jun. 1926, K. Watanabe 
(MAK); Fukuroi, 19 May 1933, K. Suzuki (KYO). — 
Aichi Pref.: Seto, 21 Jun. 1942, H. Kato 3253, (KYO, ho- 
lotype); Toyokawa, 26 Jul. 1949, K. Torii 1939 (KYO); 
Toyohashi, 22 May 1950, lnagaki (TI). —Gifu Pref.: 
Yaotsu, 18 Jul. 1931, K. Shiota (KYO). —Mie Pref.: Su- 
zuka, 7 Jul. 1936, S. Matsuyama , (KYO); Kiho, 29 May 
1999, K. Seto 50644 (OSA). —Wakayama Pref.: Nachi- 
katsuura, 18 May 1992, S. Goto , (KYO); 6 Aug. 1995, K. 
Seto 44456 { OSA); Shingu, 28 Jul. 1956, Y. Ogawa ( KYO); 
7 Oct. 1989, R. Miura (KYO); Kozagawa 21 Jul. 1976, S. 
Goto (KYO); —Hyogo Pref.: Kobe, 25 Aug. 1935, T. Ko- 
bavakawa (KYO); 27 Jul. 1937, T. Kobayakawa (KYO); 
Jul. 1937, E. Ishikawa (TI); 17 Jun. 1994, 5. Miyake 3027 
(MAK). Shikoku. Tokushima Pref.: Kainan, 10 Jun 1974, 
T. Yamanaka 66545 (TUS). —Kochi Pref.: Aki, 12 Jun. 
2007, N. lnagaki 0185122 (MBK); Konan, 2 Jun. 1973, T. 
Yamanaka 64516 (TI); 7 Jun. 1984, T. Yamazaki (TI), Su- 
saki, Jun. 1887, E. Yoshinaga (MAK); 20 May 1974, G. 
Murata et al. 310 (KYO); Yusuhara, 15 May 1973, T. 
Yamawaki M73-019 (MBK); Kubokawa, 9 Aug. 1913, H. 
Yamamoto (MAK); Mihara, 25 Jun. 2003, M. Tashiro et 
al. (MBK); Tosashimizu, 4 Jun. 1974, T. Yamanaka 66506 
(TUS). 

We are grateful to Messrs. A. Narita, M. Saji and T. 
Orime for their field assistance. Thanks are also due to 
Messrs. N. Tada and Y. Yamane for their technical sup¬ 
port. 
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Appendix 1. Mean ± standard deviation of 12 floral and three leaf characters in eight populations of V. katoi-V. yamanakae 
complex. Abbreviations of characters are given in the text. 


Population PL CLD CLL CLW CND CNL GYD GYL 

Code (mm) (mm) (mm) (mm) (mm) (mm) (mm) (mm) 


K 1 

10.3 

± 

2.6 

11.4 

± 

1.3 

4.7 

± 

0.5 

2.9 

± 

0.4 

3.0 

± 

0.3 

1.4 

± 

0.2 

1.6 

± 

0.1 

2.5 

± 

0.2 

K 3 

7.8 

± 

2.4 

9.7 

± 

1.4 

4.4 

± 

0.7 

2.1 

± 

0.3 

1.8 

± 

0.1 

0.7 

zb 

0.1 

1.3 

zb 

0.1 

1.5 

± 

0.1 

K 4 

15.7 

± 

3.1 

13.2 

± 

1.0 

5.8 

± 

0.4 

3.3 

± 

0.1 

3.1 

± 

0.2 

1.3 

zb 

0.1 

1.6 

zb 

0.1 

2.3 

zb 

0.1 

K 5 

21.4 

± 

5.3 

15.3 

± 

1.3 

6.8 

± 

0.5 

4.2 

± 

0.3 

3.3 

± 

0.2 

1.3 

zb 

0.1 

1.8 

± 

0.1 

2.6 

± 

0.2 

Y 1 

12.1 

+ 

3.4 

14.0 

± 

1.5 

7.2 

± 

1.0 

4.5 

± 

0.9 

3.9 

± 

0.3 

1.6 

± 

0.3 

1.7 

± 

0.1 

2.6 

zt 

0.2 

Y 2 

8.5 

± 

3.7 

14.8 

± 

1.3 

7.2 

zb 

0.9 

4.0 

± 

0.6 

3.9 

zb 

0.3 

1.6 

zb 

0.2 

1.6 

zb 

0.1 

2.7 

± 

0.2 

Y 3 

9.9 

± 

2.2 

13.9 

± 

1.5 

6.2 


0.7 

3.3 


0.3 

3.1 

zb 

0.3 

1.2 

zb 

0.1 

2.1 

zb 

2.7 

2.4 

± 

0.2 

Y 4 

12.7 

± 

4.5 

14.6 

± 

1.8 

7.1 

il 

1.1 

3.5 

± 

0.4 

3.6 

zb 

0.3 

1.6 

zb 

0.2 

1.8 

zb 

0.1 

2.6 

zb 

0.4 
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Appendix 1. Continued. 


Population 

Code 

COL 

(mm) 

COW 

(mm) 

POL 

(mm) 

POW 

(mm) 

LBL 

(mm) 

LBW 

(mm) 

PET 

(mm) 

K1 

0.30 ±0.03 

0.15 ±0.02 

0.26 ±0.02 

0.16 ±0.01 

14.1 ±2.8 

6.1 ± 1.2 

13.4 ±3.4 

K3 

0.23 ± 0.02 

0.12 ±0.01 

0.20 ±0.02 

0.14 ±0.01 

11.0 ± 2.8 

4.2 ± 0.9 

9.6 ±2.7 

K4 

0.25 ± 0.01 

0.14 ±0.01 

0.26 ± 0.02 

0.17 ±0.01 

11.3 ± 1.6 

4.5 ± 1.1 

10.9 ± 1.5 

K.5 

0.32 ± 0.04 

0.17 ±0.02 

0.28 ±0.03 

0.18 ±0.01 

15.4 ± 1.5 

5.9 ±0.8 

12.2 ± 4.5 

Y1 

0.28 ±0.02 

0.14 ±0.01 

0.27 ±0.03 

0.16 ±0.02 

14.0 ±2.2 

5.3 ±0.8 

15.0 ±5.0 

Y2 

0.27 ± 0.02 

0.13 ±0.01 

0.27 ±0.03 

0.17 ±0.02 

15.9 ± 1.7 

5.8 ±0.9 

15.8 ±3.0 

Y3 

0.26 ± 0.02 

0.14 ±0.02 

0.29 ± 0.03 

0.18 = 0.02 

14.7 ±2.7 

6.9 ± 1.3 

11.5 ±4.5 

Y4 

0.29 ±0.02 

0.13 ±0.01 

0.29 ± 0.02 

0.17 = 0.03 

14.0 ±2.2 

4.3 ± 0.6 

14.4 ±2.9 


Appendix 2. Allele frequencies for six loci in the Vincetoxicum katoi-vamanakae complex. 




K1 

K2 

K3 

K.4 

K5 

Y1 

Y2 

Y3 

Y4 

Aat-l 

a 

0.208 

- 

0.083 

- 

0.060 

0.120 

0.890 

0.375 

0.173 


b 

- 

0.154 

0.150 

- 

0.040 

- 

- 

0.038 

- 


c 

0.542 

0.769 

0.750 

0.700 

0.900 

0.880 

0.110 

0.588 

0.827 


d 

0.250 

0.077 

0.017 

0.300 

- 

- 

- 

- 

- 

Est-1 

a 

- 

- 

- 

0.300 

- 

0.160 

- 

0.138 

0.269 


b 

- 

- 

- 

0.700 

- 

0.840 

1.000 

0.862 

0.731 


c 

0.688 

0.615 

0.661 

- 

- 

- 

- 

- 

- 


d 

0.313 

0.385 

0.339 

- 

1.000 

- 

- 

- 

- 

Lap 

a 

0.042 

- 

0.033 

- 

0.040 

- 

- 

0.009 

- 


b 

0.604 

0.538 

0.617 

0.800 

0.960 

0.480 

0.350 

0.824 

0.385 


c 

0.146 

0.192 

0.217 

- 

- 

0.520 

0.440 

0.111 

- 


d 

0.083 

0.115 

0.133 

0.200 

- 

- 

0.210 

0.056 

0.615 


e 

0.125 

0.154 

- 

- 

- 

- 

- 

- 

- 

Mdh-1 

a 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

Pgi-1 

a 

0.188 

0.192 

- 

- 

0.880 

0.080 

- 

0.037 

- 


b 

0.563 

0.808 

0.823 

0.700 

0.120 

0.920 

1.000 

0.935 

0.981 


c 

0.146 

- 

0.081 

- 

- 

- 

- 

- 

- 


d 

0.104 

- 

0.097 

0.300 

- 

- 

- 

0.028 

0.019 

Skd 

a 

0.021 

0.154 

- 

- 

- 

- 

- 

- 

- 


b 

0.125 

- 

0.034 

0.600 

0.520 

0.320 

0.110 

0.236 

0.152 


c 

0.146 

0.231 

- 

- 

- 

0.020 

- 

- 

- 


d 

0.063 

0.308 

0.724 

0.100 

- 

0.080 

- 

0.179 

0.261 


e 

0.083 

0.308 

0.241 

0.300 

0.460 

0.580 

0.890 

0.585 

0.587 


f 

0.542 

- 

- 

- 

0.020 

- 

- 

- 

- 


g 

0.021 

- 

- 

- 

- 

- 

- 

- 

- 
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